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Abstract 



We consider an elliptic eigenvalue problem with a fast cellular flow of amplitude A, in a two- 
dimensional domain with L 2 cells. For fixed A, and L — > oo, the problem homogenizes, and has 
been well studied. Also well studied is the limit when L is fixed, and A — >• oo. In this case the 
solution equilibrates along stream lines. 

In this paper, we show that if both A — > oo and L — > oo, then a transition between the 
homogenization and averaging regimes occurs at A m i 4 . When A>i 4 , the principal Dirichlet 
eigenvalue is approximately constant. On the other hand, when A<I 4 , the principal eigenvalue 
behaves like a(A)/L 2 , where a (A) « y/AI is the effective diffusion matrix. A similar transition 
is observed for the solution of the exit time problem. The proof in the homogenization regime 
involves bounds on the second correctors. Miraculously, if the slow profile is quadratic, these 
estimates can be obtained using drift independent LP — > L°° estimates for elliptic equations with 
an incompressible drift. This provides effective sub and super-solutions for our problem. 

1 Introduction 

Consider an advection diffusion equation of the form 



where A is the non-dimensional strength of a prescribed vector field v(x). Under reasonable as- 
sumptions when A — > oo, the solution ip becomes constant on the trajectories of v. Indeed, dividing 
(11. ip by A and passing to the limit A — > oo formally shows 



which, of course, forces (p to be constant along trajectories of v. Well known "averaging" results [8] 
[T5ll20| study the slow evolution of ip(t, x) across various trajectories. 

On the other hand, if we fix A = 1, classical homogenization results [2],ll2|,f20] determine the long 
time behavior of solutions of (|l.lj) . For such results it is usually convenient to choose e < 1 small, 
and rescale (jl.ip to time scales of order 1/e 2 , and distance scales of order 1/e. This gives 
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d t ip + Av(x) • Vip - Atp = 0. 



(1.1) 



v(x) ■ Vip = 0, 




(1.2) 
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Assuming v is periodic, and that the initial condition varies slowly (i.e. <p £ (x, 0) is independent of 
e), standard homogenization results show that (p e — > <p, as e — > 0. Further, (p is the solution of the 
effective problem 

^=V-(aV<p), (1.3) 

and a is the effective diffusion matrix, which can be computed as follows. Define the correctors xi> 
. . . , Xn to be the mean-zero periodic solutions of 

- Axj +v(x) -Vxj = -Vj(x), j = l,...,n. (1.4) 

Then 



+ ^Xi-^Xjdx, i,j = l,...,n. (1.5) 
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Q is the period cell of the flow v(x), and 5{j is the Kronecker delta function. 

The main focus of this paper is to study a transition between the two well known regimes 
described above. To this end, rescale (II. lh by choosing time scales of the order 1/e 2 and length 
scales of order 1/e. This gives 

9t<Pz,A-\ V ( - ) • Vp e ,A - A(fe,A = 0, (1.6) 



e ve 

where A 3> 1 and e <C 1 are two independent parameters. Of course, if we keep e fixed, and send 
A — > oo, the well known averaging results apply. Alternately, if we keep A fixed and send e — > 0, 
we are in the regime of standard homogenization results. The present paper considers (|1.6p with 
both e — > and A — > oo. Our main result shows that if v is a 2D cellular flow, then we see a sharp 
transition between the homogenization and averaging regimes at A ~ 1/e 4 . 

Before stating our precise results (Theorems 11.11 and 1 1 . 2 1 b elow) . we provide a brief explanation as 
to why one expects the transition to occur at A ~ 1/e 4 . For simplicity and concreteness, we choose 
the stream function H(x\,X2) = ^ sin(7rxi) sin(7TX2), and define v{x\,X2) = (— d-iH, d\H). Even in 
this simple setting, to the best of our knowledge, the transition from averaging to homogenization 
has not been studied before. 

First, for any fixed A, we let a (A) = ((jy(A)) denote the effective diffusion matrix obtained in 
the limit e — > (see |18| for a comprehensive review). If xf is the mean zero, 2-periodic solution to 

- A X f + Av(x) ■ V Xj = -Avj(x), for j E {1, 2}, (1.7) 

then the effective diffusivity (as a function of A) is given by fjl .5j) . As A — > oo, the behaviour of 
the correctors xf is weu understood [5l[6l[ini[IZl[l9l[2Tl[2ll[25] . Except on a boundary layer of order 
1/y/A, each of the functions Xji x ) + x j become constant in cell interiors. Using this one can show 
(see for instance [HE]) that asymptotically, as A — > oo, the effective diffusion matrix behaves like 

a(A) = a VAI + o(VA) (1.8) 

Here I is the identity matrix, and o~q > is an explicitly computable constant. Consequently, if we 
consider (|1.6j) . with the Dirichlet boundary conditions on the unit square, we expect 

(p £) A.(x, t) exp(— ao\^At), as t — > oo, for small e. (1-9) 

On the other hand, if we keep e fixed and send A — > oo, we know [8, 25 that ip becomes constant 
on stream lines of H. In particular, because of the Dirichlet boundary condition on the outside 
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boundary, we must also have <p = on the boundary of all interior cells. Since these cells have side 
length s, we expect 

(p Et A{x,t) ~ exp(— ir 2 t/e 2 ), as t — > oo for large A. (1-10) 

Matching (|1.9p and (jl.lOp leads us to believe \f~A ~ 1/e 2 marks the transition between the two 
regimes. 

With this explanation, we state our main results. Our first two results study the averaging 
to homogenization transition for the principal Dirichlet eigenvalue. Let L be an even integer, 
D = [-L/2, L/2] x [— L/2, L/2] be a square of side length L, and A > be given. We study 
the principal eigenvalue problem on D 

-A<p + Av ■ Vy? = X(p in D 

ip = ondD (1.11) 
if > in D, 

as both L, A — > oo. We observe two distinct behaviors of A with a sharp transition. If A L 4 , then 
the principal eigenvalue stays bounded, and can be read off using the variational principle in [3] in 
the limit A — > oo. This is the averaging regime, and exactly explains (jl.lOj) . On the other hand, if 
A <C L 4 , then the principal eigenvalue is of the order a(A)/L 2 . This is the homogenization regime, 
and when rescaled to a domain of size 1, exactly explains (|1.9p . Our precise results are stated below. 



Theorem 1.1 (The averaging regime). Let ip = <pl,A ^ e the solution of (jl.lip and A = \l A be the 

principal eigenvalue. If A — >■ oo, and L = L(A) varies such that 



lim inf r9l " ^ > 0, (1.12) 

A^oo L 2 log A log L v ; 

then there exist two constants Ao, Ai, independent of L and A, such that 

< A < Xl,a ^ Ai < oo (1.13) 

for all A sufficiently large. 

Theorem 1.2 (The homogenization regime). As with Theorem let ip = <pl,a be the solution 
of (II. lip and A = Xl,a be the principal eigenvalue. If L —> oo, and A = A(L) varies such that 

-L 4 ~ a ^A^ cL 4 ~ a , for some a > 0, (1.14) 
c 

then there exists a constant C = C(a,c) > 0, independent of L and A, such that 

~§ * ^a < C§- (1.15) 

for all L sufficiently large. 

Remark. In the special case where A = L@, for f} > 4, assumption (|1 . 12[) is satisfied, and consequently 
the principal eigenvalue remains bounded and non-zero. For (3 < 4, assumption (|1.14|) is satisfied 
and the principal eigenvalue behaves like that of the homogenized equation. 
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In the averaging regime (Theorem II .ip . the proof of the upper bound in f) 1 . 13[) follows directly 
using ideas of [3]. The lower bound, however, is much more intricate. The main idea is to control 
the oscillation of ip between neighbouring cells, and use this to show that the effect of the cold 
boundary propagates inward along separatrices, all the way to the center cell. The techniques used 
are similar to [3[T6]. The main new (and non-trivial) difficulty in our situation is that the number 
of cells also increases with the amplitude. This requires us to estimate the oscillation of <p between 
cells in terms of energies localised to each cell (Proposition 12. 4| below). Here the assumption that 
L is not too large comes into play. Finally, the key idea in the proof is to use a min-max argument 
(Lemma 12. 5} below) to show that <p is small on the boundaries of all cells. 

Moreover, once smallness on separatrices is established, our proof may be modified to show that 
under a stronger assumption 

-/A 

(1.16) 



lim inf — ^ 

A->oo L A log A log L 



-oo, 



we have a precise asymptotics 



lim \l A = inf 



iViol 



w e Hl(Q), 



w 



1, and w ■ Vv = 



Q 



(1.17) 



where Q is a single cell. This is the same as the variational principle in [3]. We remark however 
that [3] only gives (|1.17p for fixed L as A — > oo. 

Turning to the homogenization regime (Theorem ll.2p , we remark first that homogenization of 
eigenvalues has not been as extensively studied as other homogenization problems. This is possibly 
because eigenvalues involve the infinite time horizon. We refer to [Hll3 l fT4" l l22 p 23] that all study self- 
adjoint problems for some results on the homogenization of the eigenvalues in oscillatory periodic 
media. The extra difficulties in the present paper come both from two sources. First, since the 
problem is not self adjoint, a variational principle for the eigenvalue is not available. Second, as A 
and L tend to oo, we don't have suitable aprori bounds because either the domain is not compact, 
or the effective diffusivity is unbounded. 

Our proof uses a multi-scale expansion to construct appropriate sub and super solutions. When 
A is fixed, it usually suffices to consider a multi-scale expansion to the first corrector. However, in 
our situation, this is not enough, and we are forced to consider a multi-scale expansion up to the 
second corrector. 

Of course an asymptotic profile, and explicit bounds are readily available [6] for the first corrector. 
However, to the best of our knowledge, bounds on the second corrector as A — > oo have not been 
studied. There are two main problems to obtaining these bounds. The first problem is appearance 
of that terms involving the slow gradient of the second corrector multiplied by A. In general, we 
have no way of bounding these terms. Luckily, if we choose our slow profile to be quadratic, then 
these terms idnetically vanish and present no problem at all! 

The second problem with obtaining bounds on the second corrector is that it satisfies an equation 
where the first order terms depend on A. So one would expect the bounds to also depend on A, 
which would be catastrophic in our situation. However, for elliptic equations with a divergence free 
drift, we have apriori L p — > L°° estimates which are independent of the drift [HIT]. This, combined 
with an explicit knowledge of the first corrector, allows us to obtain bounds on the second corrector 
that decay when A <C L 4 . 

The sub and super solutions we construct for eigenvalue problem are done through the expected 
exit time. Since these are interesting in their own right, we describe them below. Let r = tl,a be 
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the solution of 



-At + Av ■ Vr = 1 in D 
t = on dD, 



(1.18) 



where v and D are as in (jl.lip . Though we do not use any probabilistic arguments in this paper, it 
is useful to point out the connection between r and diffusions. Let X be the diffusion 



dX t = -Av(X t ) dt + V2dW t 



(1.19) 



where W is a standard 2-dimensional Brownian motion. It is well known that r is the expected exit 
time of the diffusion X from the domain D. Numerical simulations of three realizations of X are 
shown in Figured) Note that for "small" amplitude (A = L 3 ), trajectories of X behave similarly to 
those of the Brownian motion. For a "large" amplitude (A = £ 4 ' 5 ), trajectories of X tend to move 
ballistically along the skeleton of the separatrices. 






o 




'fj' ! 3' 1 


o 


o 








© 






o 




Q 







: 






Q 




Q 





(a) Small amplitude (A — L 3 ) (b) Large amplitude (A = L 4 ' 5 ) 

Figure 1: Trajectories of three realizations of the diffusion (11.19j) , 




(a) Small amplitude (A = L 3 ) (b) Large amplitude (A — L J ) 



Figure 2: A contour plot of r(x,y). 

Similar to the eigenvalue problem, the behaviour of r is described by two distinct regimes with 
a sharp transition. If A 3> L 4 , then the stirring is strong enough to force the diffusion X to exit 
D almost immediately along separatrices. In this case, we show that r — > on separatrices, and 
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is bounded everywhere else above by a constant independent of A and L. On the other hand, if 
A <C L 4 , then the stirring is not strong enough for the effect of the cold boundary to be felt in 
the interior. In this case, it takes the diffusion X a very long time to exit from D, and r — > oo as 
A,L—> oo. A numerical simulation showing r in each of these regimes is shown in Figure EJ The 
precise results £11*6 clS follows. 



Theorem 1.3 (The averaging regime). Let r = tl,a be the solution to (11.180 . Let A — > oo, and 
suppose L = L(A) varies such that (|1.12p is satisfied. There exists a constant C, independent of A, 
L, such that for all A sufficiently large 

L 2 

t(x) 2 ^ C—= log A log L, whenever H(x) = 0. 
v A 



Consequently, if H(x) = 0, then t(x) — > as A — )• oo, and \\t\\ Loo ^ is bounded uniformly in A. 

Theorem 1.4 (The homogenization regime). As with Theorem \1.3l let r = tl,a be the solution 
of (I1.18P on the square D = [-L/2, L/2] x [— L/2, L/2]. Suppose now L — > oo, and A = A(L) varies 
such that ()1.14p is satisfied, for some fixed a 6 (0,4). Then, for any 5 > 0, there exists a constant 
C = C(5,a,c) > 0, independent of A, L, such that 



L 



L 



L 



C —= ^ t{x) ^ C—=, whenever \x\ ^ (1 — 5) — 



(1.20) 



for all L sufficiently large. Consequently, t — > oo as L — >■ oo, uniformly on compact sets. 



The proof of Theorem 11.31 is similar steps to that of Theorem 11.11 For the proof of Theorem 11.41 
as mentioned earlier, we need to perform a multi-scale expansion up to two correctors, and choose 
the slow profile to be quadratic. When the domain is a disk, a quadratic function is exactly the 
solution to the homogenized problem! This gives us a sharper estimate for r. 



Proposition 1.5. Let Bl be a disk of radius L, and ta,l be the solution of (|1.18p in B^. If A and 

L satisfy the assumptions in Theorem \l-4\ then 



t(x) 



1 



2tr(a(A)) 



L 1 



^ c 



L 



A 1 / 4 



(1.21) 



where c > is independent of A and L. Here <j(A) is the effective diffusion matrix and tr(a(A)) 
denotes the trace of this matrix. 



Remark 1.6. Note that right hand side of f[ 1 . 2 1 1) tends to infinity as A, L — > oo. However, by (jl.8p the 
terms on the left are of order L 2 / \/~A, which dominates the right hand side. Thus (ll.2ip immediately 
implies (fL20|) . 

By fitting a disk inside, and outside a square, Proposition 11.51 quickly implies Theorem 11.41 
Further, since it is well known that the principal eigenvalue is bounded below by the maximum 
expected exit time, the lower bound in Theorem 11.21 also quickly follows from Proposition 11.51 The 
upper bound is a little more technical, however, also uses Proposition 11.51 as the main idea. 

We mention that we have chosen to use the particular form of the stream-function H{x\,X2) = 
— sin(-7rxi) sin(-7TX2) simply for the sake of convenience. All our results may be generalized to other 
periodic flows with a cellular structure without any difficulty. We also believe that for other flows the 
transition from the averaging to the homogenization regime happens when the effective diffusivity 
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(7 (A) balances with the domain size L. That is, when the "homogenized eigenvalue" a(A)/L 2 is of 
the same order as the "strong flow" eigenvalue: 

lim „ 1. (1.22) 

We leave this question for a future study. 

This paper is organized as follows. The averaging regime is considered in Sections [2] and [3j The 
former contains the proof of Theorem 11.11 and the latter of Theorem 11.31 The rest of the paper 
addresses the homogenization regime. The key step here is Proposition 11.51 proved in Section HJ 
From this, Theorem 11.41 quickly follows, and the proof is presented the same section. Theorem 11.21 
is proved in Section [H 
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2 The eigenvalue in the strong flow regime 

In this section we present the proof of Theorem 11.11 First, we discuss the proof of the upper 
bound in (|1.13p . followed by the proof of the corresponding lower bound, and, finally, of the limiting 
behavior in (j!.17jl . 

2.1 The upper bound 

The upper bound for A in (|1.13p follows directly from the techniques of |3j. We carry out 
the details below. Following [3], given any test function w G Hq(D), and a number a > 0, we 
multiply (jl.lip by w 2 /(f + a) and integrate over D to obtain 

. f w 2 ld f w 2 Aip , f w 2 „ , 

A / 1-=- / -+A / v-Vip. (2.1) 

J D tp + a J D tp + a J D tp + a 

For the first term on the right, we have 

w 2 Atp f ^ /2w((p + a)Vu) — w 2 V(p 

D ip + a J D V (f + a ) 2 

2 f \wVip- (ip + a)Vw\ 2 f 2 
Vw\ - / — ■ — -2 < / \Vw\ . 



For the second term on the right of (J2JJ) we have, since u is incompressible, 

/ v-Vtp= / w 2 v ■ V ln(tf + a) = — 2 / ln(<^ + a)w(v ■ Vw). 

JD<P + a ' Jd ' Jd 



Hence, equation (|2.ip reduces to 

,,2 



A 



f ^ f \V W \ 2 -2A [ ln(ip + a)w{vVw). (2.2) 

J D ip + a J D J D 
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Now, choose w to be any Hq(D) first integral of v (that is, v ■ Vro = 0). Then, equation (|2,2p 
reduces to 

A f ^ < [ |VH 2 . 
J D if + a J D 

Upon sending a — > 0, the Monotone Convergence Theorem shows 

A / w 2 < / | V^| 2 

JD JD 

for any Hq(D) first integral of u. Choosing w = H(x), which, of course, does not depend on L, we 
immediately see that 

A^ ( [ H 2 ) 1 [\VH\ 2 = (l 2 [ H 2 ] 1 L 2 I \VH\ 2 =( I H 2 \ % f \VH\ 2 , 
\Jd ) JD V JQo J JQo \JQo / JQo 

where Qq is any cell in D. This gives a finite upper bound for A that is independent of L and A. 

2.2 The lower bound 

The outline of the proof is as follows. The basic idea is that if the domain size L is not too 
large, and the flow is sufficiently strong, the eigenfunction ip should be small not only near the 
boundary dD but also on the whole skeleton of separatrices inside D. Therefore, the Dirichlet 
eigenvalue problem for the whole domain D is essentially equivalent to a one-cell Dirichlet problem, 
which gives the correct asymptotics for the eigenvalue for A large. To this end, we first estimate the 
oscillation of ip along a streamline of v inside one cell that is sufficiently close to the separatrix, and 
show that this oscillation is small: see Lemma l2.2i Next, we show that the difference of the values 
of ip on two streamlines of v (sufficiently close to the separatrix) in two neighbouring cells must be 
small, as in Lemma 12.31 below. These two steps are very similar to those in [7], and their proofs are 
only sketched. 

Now, considering the 'worst case scenario' of the above oscillation estimates, we obtain a point- 
wise upper bound on <p on streamlines of v near separatrices in terms of the principal eigenvalue 
A, and ||</?||^,2: see Lemma 12.51 Next, we show that the streamlines above enclose a large enough 
region to encompass most of the mass of ip 2 . Finally, we use the drift independent apriori estimates 
in [HE] to obtain the desired lower bound on A. 

A streamline oscillation estimate 

The basic reason behind the fact that the eigenfunction is constant on streamlines is the following 
estimate, originally due to S. Heinze flO] . 

Lemma 2.1. There exists a constant C > so that we have 

J \v ■ V^\ 2 dx ^jj \V^\ 2 dx = ^|M|| 2 . (2.3) 

Proof. Let us use the normalization ||</?||l 2 = 1- We multiply (jl.lip by v ■ and integrate over 
Qi. This gives 

A / \v\7tp\ 2 = X / (p{v-Vip)+ / Atp(vVtp). 

JD JD JD 

Notice that 

p(v-V V ) = ^J D VV{<p 2 ) =0, 
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since v ■ v = on dQi and V • v = 0. Similarly, we have, as v ■ Vy? = on dD: 

2 2 

/ Aip(v ■ V99) = - V / djip(djv -V<p) -^2 dj<p(v -Vdj<p)+ \ {v ■ V<p) (v • V<p) dS 

JD ~[ JD ~[ J D J 3D 

^l|V W |U (Qi) ||V^||| 2(D) -^ f V-(v(d j pf)=C\ 

3=1 JD 

and consequently we obtain (|2.3p . □ 

The next lemma bounds locally the oscillation on streamlines in terms of the L 2 -norm of v • V92. 

Lemma 2.2. Let Qi be any cell. For any 5o > 0, there exists Tj C Qi such that Ti is a level set of 
H, \H(Ti)\ E (5 , 2S Q ), and 

sup \ip{x x ) - <p{x 2 )\ 2 ^C^log(j-) [ \vVtp\ 2 . (2.4) 

for some constant C independent of A,L,5q. 

We will see that 5q is the 'width' of the boundary layer, and will eventually be chosen to be 
5 « 1/VA. 

Proof. The proof is straightforward and similar bounds have already appeared in [71ll6j. We sketch 
the details here for convenience. First we introduce curvilinear coordinates in the cell Qi. For this, 
let (xi,yi) be the center of Qi, and @i be the solution of 



V@i-VH = in Qi-{(xi + t, yi ) I i>0} 

Qi(x, y) = tan -1 ( — — — ) on dQi. 

\X-XiJ 

As usual, we extend to Qi by defining it to be (or 2tt) on {(xj + 1, yi) \ t ^ 0}. 
In the coordinates (h, 8) given by the functions H and 0j, it is easy to check that 

dp v ■ Vp 



(2.5) 



de |v©i||vjy| " 

Assume, for simplicity, that H ^ on Qi. Then for any h E (8q, 25q), we have 

/ \ 2 

de \ 



sup|^(Mi)-<KM 2 )| 2 < / \v-V<p\ 
61,62 \J{H=h} |VHi||V-H| 

,2 d6 f fit! 



J {H =h} 1 |ve<||v#| J{H=h} |ve,||VF 

ZCln(j-) [ \v-Vp\ 2 d9 



SoJ J{H=h} |ve,||VF| 



The last inequality follows from the fact that 



(«.;,) |ve,||VH 



'"' <£7tai. 
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as the length element along the contour dl = d9/\VQ\. Now integrating over (5q,26q), we get 

/•2<$o i r 

I sup\ip(h,6 1 )-ip(h,6 2 )\ 2 dh^Cln— \v ■ Vp\ 2 
J 5 0i,02 ^0 Jq, 

and (12.41) follows from the mean value theorem. 



□ 



Variation between neighboring cells 

Now, we consider two streamlines on which the solution is nearly constant and estimate the 
possible jump in the value of tp between them. 

Lemma 2.3. Let Qi and Qj be two neighbouring cells, Tj C Qi, Tj C Qj the respective level sets 
from Lemma \2. 6 A and let hi = H(Ti), hj = H(Tj). Then there exists Xi S Tj and Xj G Tj such that 



\<p(xi) - (p(xj)\ 2 ^ C5 



! |Vvf + cllog(-M / \v-Vp\ 2 . 



(2.6) 



Proof. Assume again for simplicity that H ^ on Qi, and Qi is to the left of Qj. Then, using the 
local curvilinear coordinates (h,6) around the common boundary between the cells Qi and Qj, as 
in (|2.5p . we have 

<p(h i ,e)-<p(h j ,o) = ^dh. 

Now, let 5\ £ (0, ^) be fixed. In the region \h\ ^ hi, and \8\ ^ 6\, we know that |V.ff| ~ 1 and 
|VO| ~ 1. Hence, we have 



9 / ,<5 i f h J 

\Vp\ 2 > C / / 



dip 



dh 



Oq |0|<5i 



However, Lemma 12.21 shows that 



sup 0)| 3 < inf |^,0)-^,0)| 2 + cilogf-i') / 



\v ■ Vp\' 



This concludes the proof of (|2.6p . 



(2.7) 
□ 



Variation between two far away cells 

For each cell Qi, we set 

at = [ (\Vip\ 2 + A\v Vp\ 2 ) dx. 
JQi V ' 

Choosing 5q = 1 /VA, Lemmas 12.21 and 12.31 immediately give the following oscillation estimate. 

Proposition 2.4. If Qi and Qj are any two cells, and Tj C Qi, Tj C Qj the respective level sets 
from Lemma UTB, then 



sup ip{xi) — inf ip{xj" 



^ C A l/4 LV^' 



line 



where the sum is taken over any path of cells that connects Qi and Qj, consisting of only horizontal 
and vertical line segments. 
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Lemma 12 . 1 1 implies that 

J>^CA|M|| 2(fl) , (2.8) 

3 

with the summation taken over all cells in D. Now, the key to the proof of the lower bound in 
Theorem 11.11 is to obtain an estimate on |Mlz,oo( r .) in terms of (Yli^) 1 ^ 2 - A direct application of 
Cauchy-Schwartz to Proposition 12.41 is wasteful and does not yield a good enough estimate. What 
is required is a more careful estimate of the 'worst case scenario' for the values of ol{. This is the 
content of our next Lemma. 

Lemma 2.5. On any cell Qi, we have 

i , v,2 ~ log 4 log Z ^ log^logL 2 

sup|y>(x)| < C A1/2 ^ ai ^C A±/2 \\\tp\\ L 2 iD) 

1 all cells 



where Ti C Qi is the level set from Lemma \2.2l 

Proof. For notational convenience, in this proof only, we will assume that D = (-L — ^, L + is 
the square of side length 2L + 1 centered (0,0), and Qi,j = {(x,y) \ x G [i — ^,i + ^),y^ [j— oji+o)) 
is the cell with center Note that in the present proof we label the cells, (and contours 

inside the cell Qij we use below) by two indices that correspond to the coordinates of the center of 
the cell. 

Let &i j denote the set of all paths of cells that join the boundary dD to the cell Qi ,j using 
only horizontal and vertical line segments. Let /3 = (/3ij) G M^ 2L+1 ^ , /3jj ^ 0, be a collection of 
non-negative numbers assigned to each cell, and denote 

QioJoW) = min Y] \fKr ( 2 - 9 ) 
We first claim there exists an explicitly computable constant C, independent of L,/3,io, jo such that 

L 

q l0 , j0 W) 2 ^ClogL (2-10) 

i,j=-L 

which is an obvious improvement over the Cauchy-Schwartz estimate applied blindly to (|2.9p . This 
improvement comes because we are taking the minimum over all such paths in (12. 9h . 
To prove (|2.10p . we define 

1 bojoV e©; OiJO (i,i)e fl 

where ®', io ^ is any collection of (possibly repeated) paths in &(i j )- Since the minimum of a 
collection of numbers is not bigger than the average of any subset, we certainly have 

W/3K<^/3,<5' loJO ) 

for any collection <5' iQ JQ . The idea is to choose such a sub-collection in a convenient way. 

We prove the claim for (io, jo) = (0,0). We choose <5' 00 to consist of (L + 1)! paths, with the 
following property. All paths stay in the upper-right quadrant. The last cell visited by all paths is 
(0,0). The second to last cell visited by (L + l)!/2 paths (half of the collection) is (1,0), and the 
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second to last cell visited by the remaining paths is (0, 1). Amongst the paths who's second to last 
cell is (1,0), we choose O o so that two thirds of these paths have (2,0) as the third to last cell, 
and one third have (1,1) as the third to last cell. Symmetrically, we choose (5 o so that amongst 
all the paths who's second to last cell is (0, 1), two thirds of these paths have (0, 2) as the third last 
cell, and one third have (1, 1) as the third to last cell. Consequently exactly {L + l)!/3 paths have 
third to last cell (2,0), exactly (L + l)!/3 paths in 0q O have third to last cell (1,1), and exactly 
(L + 1)1/3 paths in & 00 have third to last cell (0, 2). 

Continuing similarly, we see that <25q can be chosen so that for any cell with i + j ^ L, 
exactly (L + l)!/(i + j + 1) paths visit the cell as the (i + j + l) th to last cell. Finally, we 

assume that all paths in <5' start on the top boundary and proceed directly vertically downward 
until they hit a cell of the form (k, L — k). 

Let us count how many times each term y / Aj appears in the averaged sum gg V g (/3, ©o,q)- Clearly, 

if i+j ^ L, then the cell appears in exactly ^"j^ paths in & . On the other hand, if i+j > L, 
then the cell appears exactly j^n^ paths. Consequently, we have 

Co(A^o) = E^Ev / A^+ r ^ T E E v^v (2.H) 

fc=0 i=0 i=0 j=L-i 

We now maximize the sum in f)2.11j) with the constraint 

L 

E &i = °- ( 2 - 1 2 ) 

Let S denote the right side of ()2.1ip . then at the maximizer of S we have 

dS 1 



and 

OS 1 



The Euler-Lagrange equations now imply that 



for i + j > L, 



, for i + j ^ L. 



Pij = 4(L + 1)2 fori + -? >L ' 



and 

Pij = ,/- x 7 x1 w for 0^ i + j ^L. 
Here 7 is the Lagrange multiplier that can be computed from the constraint (|2.12|) : 



L-l 

l) 2 7 ^4(7 + l) 2 



7 (L + l)(L + 2) , ^ (j + 1) 

a. 



4(L + 1) 2 '^4(j 

It follows that 

7 = <ri{L), 
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with 7(L) = 0(1/ log L) as L — > +00. Hence, for the maximizer we get 

and thus (|2.10p holds for io,jo = (0,0). However, it is immediate to see that the previous argument 
can be applied to any cell considering appropriate collection of paths that say up and to the right 
of (io,jo), whence (|2.10p holds for all («o> Jo)- 

With (|2.1U|) in hand, we observe that Proposition 12.41 implies 

11 1,2 ^ n l ogA 2 log A A A logAlogL 2 

^Loo(r ioJo ) < c^jHojoW < c-^iogL 2^ 2^ «m < c — — A IMIl2> 

where the last inequality follows from Lemma 12.11 This concludes the proof. □ 

Our next step shows that the mass of ip 2 in the regions enclosed by the level sets Tj is comparable 
to IMl£a(ij). 

Lemma 2.6. Let Qi be a cell, and Tj C Qi the level set from Lemma \2.SX Let hi = H(Ti), and 
Si = Qi H {\H\ < \hi\} be a neighbourhood of dQi. Let Q[ = Qi — Si. Then, for A sufficiently large, 
we have 



Proof. For any cell Qi, the Sobolev restriction theorem shows 

j H _ \>p(h, e)\ 2 — — = IMl£a (H -i (h)nQ4) < c\\^f H t m = c (on + 



L 2 (Qi)l ■ 



Thus, using curvilinear coordinates with respect to the cell Qi, and assuming, for simplicity, that 
hi = H(Ti) > 0, gives 



hi r2ir 



^hs l ) = I „ / y^ e ) 2 \x?Pi\\x?m dedh 



h=o Je=o 



h=o \fh 
Summing over all cells gives 

£lM&(*} < C ^ (WV^WhiD) + Ml*iB)) < Cy% (1 + A) 11 ^ Cv^ll V\\l 2 {D) 



where the last inequality follows using the upper bound in 13[) which was proved in Section 12.11 
Since 5o — > as A — > 00, and 

IMIl2(d) = ElMli^) + £MIl 2 (qo' 

i i 

inequality (I2.13P follows. □ 
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Our final ingredient is a drift independent LP — > L°° estimate in [I]. We recall it here for 
convenience. 

Lemma 2.7 (Lemma 1.3 in [2]). Let Q C K rf be a domain, w be divergence free, and 9 be the solution 
to 

f -A9 + W-V6 = f inil, 
| 8 = ondtt, 

with f 6 L p (Vt) for some p > d. There exists a constant c = c{VL,d,p) > 0, independent of w, such 
that \\9\\ Loo ^ c||/|| iP . 



We are now ready to prove the lower bound in Theorem 11.11 

Proof of the lower bound in Theorem 11.11 Using the notation from Lemma 12. 6| define D' = 
\J i Q\, and let Qj be a cell such that |MIlo°(q' ) = IMI l°°(D')- Then 

2 ^ — ^ 2 2 2 

IMIl 2 (D') = 2JMIl 2 (QJ) ^ L W^L^iQ'^ ( 2 - 14 ) 
i 

and it follows from Lemmas 12.51 and 12.71 that 

IMIl°°(q:) ^ C (a|MIl°°(Q^) + IMIz,°°(I\,-)) 

< C U\\<p\\ L ~wj + -^(logAlogL) 1 / 2 ^^!!^^)) 

*s c (mMl^ + (log AiogL^Vxyw^,^ 

where the last inequality follows from Lemma 12.61 Consequently, 

1 , IMli-(Q') A 1 ' 2 A 1 ' 2 
A > — - or \> ~ —z-TT--. r-, =- ^ 



2C \\(p\\h(D>) 2Cl °sAlogL ' 2CL 2 logAlogL 

where the last inequality follows from equation (|2.14p . This proves the lower bound on A in Theo- 
rem [TTTJ □ 

3 The exit time in the strong flow regime 

In this section we sketch the proof of Theorem 11.31 The techniques in Section 12.21 readily show 
that oscillation of r on stream lines of v becomes small. Now, the key observation in the proof of 
Theorem 11.31 is an explicit, drift independent upper bound on the exit time. We state this below. 

Lemma 3.1 (Theorem 1.2 in |llj). Let Vt C M. d be a bounded, piecewise C 1 domain, and u : Q. — > W 1 

a C 1 divergence free vector field tangential to dQ,. Let t' be the solution to 

[ - At' + u • Vr' = 1 in CI 

i , (3-1) 

\ t' = ondQ, 

Then for any p € [1, oo], 

II — ' II < || 

ll T \\Lp(Q) ^ ll T rlli>(B)' 

where B C W 1 is a ball with the same Lebesgue measure as Q., and r' r is the (radial, explicitly 
computable) solution to (13. ip on B with u = 0. 
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With this, we present the proof of Theorem 11.31 

Proof of Theorem 11.31 Following the same method as that in Section 12.21 we obtain (analogous 
to Lemma l2.5p 

i / m2 ^ Clog^logL ^ 
sup|r(x)| — > a h (3.2) 

^ VA all cells 

with ctj now equal a« = Jq.|Vt| 2 , where Qi is the i th cell. The sets Tj C Qi appearing in (|3.2p are 
level sets of H on which the oscillation of r is small (analogous to Lemma l2.2h . 
Now observe that 

E <* = /J Vr 

all cells 

and so (13.21) reduces to 



2 

T, 

D J D 



i t \\2 ^ ClogAlogL . 
sup|r(x)| ^ — / r. (3.3) 



Letting Q\ be the region enclosed by Tj, we obtain (similar to Lemma |2. 6 



/ T < 2 E /, 



(3.4) 

all cells " ^ i 



for large enough A. By Lemma 13.11 we see 

/ T^('[\ + || T" 

and hence 



L 



1/2 



' D ^ ' \JD 

Solving the above inequality quickly yields 

r sC c( L 2 + -^LlogAlogL) ^ CL 2 
d V VA 



where the second inequality above follows from the assumption (|1.12j) . Substituting this in (|3.3p 
immediately shows that 

L 2 

\\ T \\ 2 L°°(ri) < C-^logyllogL. 

Now, to conclude the proof, we appeal to Lemma I3T1 again. Let S = D — Li iQ[ be the (fattened) 
skeleton of the separatrices. Observe that \S\ ^ which, by assumption (|1.12p . remains bounded 
uniformly in A. Consequently, by Lemma 13. 11 

L 2 

||r|| ioo(5) < C\S\ 2 + \\r\\ Loo{ds) < C-^=log^logL, 



which immediately yields the desired result. □ 
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4 The exit time in the homogenization regime 

Exit time from a disk. 

The key step in our analysis in the homogenization regime is Proposition [T3J and we begin with 
it's proof. The idea of the proof is to construct good sub and super solutions for the exit time 
problem in a disk of radius one. Let r be the solution of (|1.18p in a ball of radius L. Let B\ be a 
ball of radius 1, and let T\(x) = t(Lx)/L 2 . Then t\ is a solution of the PDE 



J -An + ALv(Lx) • Vri = 1 in B%, 
| t\ = on dB\. 

We begin by constructing an approximate solution fi, by defining 



(4.1) 



fi(x) = tio(x) + -m(x,y) + jpTi2(y), (4.2) 



where y = Lx is the 'fast variable'. We define t\q explicitly by 

1 - Ixl 2 

tw(x) = ±- L , (4.3) 



and obtain equations for m and t\i using the standard periodic homogenization multi-scale expan- 
sion. Using the identities 

V = V x + LV y and A = A^ + 2LV X ■ V y + L 2 A y 

we compute 

-Af + ALv • Vf = - A^io + ALv ■ V x t 10 

1 



+ -( - A x t u - 2LV X ■ V y T U - L 2 A y n 



+ LAv ■ V x rii + L 2 Av ■ V y Tn) 



+ (-L 2 A y T 12 + L 2 Av ■ V y T 12 ) . 

We choose rn to formally balance the O(L) terms. That is, we define m to be the mean-zero, 
periodic function such that 

- A y Tu + Av ■ Vym = -Av(y) ■ V x t w . (4.4) 

We clarify that when dealing with functions of the fast variable, we say that a function 9 is periodic 
if 9(yi + 2, 2/2) = 6 1 , 2/2 + 2) = 9(yi,y 2 ) for all (yi, y 2 ) £ M 2 . This is because our drift v is periodic, 
with period 2 in the fast variable, and each cell is a square of side length 2, in the fast variable. 

Now we choose t\ 2 to formally balance the 0(1) terms. Define t\ 2 to be the mean-zero, periodic 
function such that 

- A y T 12 + Av ■ V y T 12 = 2V X ■ VyT U - A (v ■ V x T X i - {v ■ V x t u )) , (4.5) 

where (•) denotes the mean with respect to the fast variable y. Observe that we had to introduce the 
term A{v ■ V x th) above to ensure that the right hand side is mean zero, to satisfy the compatibility 
condition. 
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We write 

m(x,y) = xi{y)d xl Ti {x) + x2(y)d X2 T 10 (x) = -xi{y)x\ -X2(y)x 2 , (4.6) 

where Xj = Xj{y)i J = 1,2 are the mean zero, periodic solutions to 

- AyXj + Av • V y xj = -Avj. (4.7) 

Using this expression for m and (I4.3P we simplify (I4.5P to 

- A. y ri2 + Av • V y T 12 = -2d yi xi - 2d V2 X2 + A(vixi + ^2X2 - (viXi) ~ (^2X2))- (4.8) 

The key observation is that with our choice of no, the right side of (|4.8p is independent of the 
slow variable. Our aim is to show that f satisfies the estimates (14. 9p and (]4.10p below. 

Lemma 4.1. There exists a positive constant cq = cq(o) independent of A, L, such that for f 
defined by (j4.2j) we have 

\f{x) - no (a) I < c L- a/4 for x G #i (4.9) 



and 

- Af + ALu(Lx) • V x f = tr(a(A)). 
iiere <r(^4) is the effective diffusion matrix, given by (jl.5p . 

We first use the Lemma to finish the proof of Proposition 11.51 



(4.10) 



Proof of Proposition 11.51 The key observation we obtain from Lemma 14. II is that, except for the 
boundary condition, the function f'(x) = f(x)/ ti(a(A)) satisfies exactly (j4. 1[) . This is a miracle 
that happens only when the domain is a disk. Then we get sub- and super-solutions for t\[x) by 
setting 



and 



t(x) 
t(x) 



1 



tr(a(A)) 
1 



f(x) + 



TlX) 



2c 

L a/A 

2c 

L a/4 



Lemma |4. II implies that 

-Ar(x) + ALv(Lx) • Vr(x) = -Ar(z) + ALv(Lx) • Vt(x) = 1. 

Further, since no 0*0 = on dB\, equation (|4.9p implies that f{x) > and r(x) < on dB\. 
Consequently, r is a super solution, and r is a sub solution of (14. ip . and hence 



1 



f(x) 



2c 

L a/4 



tv(a(A)) 

Rescaling to the ball of radius L, we see 

L 2 



s$ n{x) < 



tr(*(A)) 



t(x) 



tv{a(A)) Tl ° \L 



L 2 



f{x) + 



4c 



2c 

L a/4 



(4.11) 



tr(a(A)) L a / 4 ' 



Now using (fL8j) and (fTTIj) we obtain (OB . 
It remains to prove Lemma 14.11 



□ 
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Proof of Lemma 14.11 By our definition of no, th, t±2, we have 



-Af + ALv ■ Vf = -A x t 10 - -A xTn - A[(v lX i) + (V2X2}} = 2 + (|Vxi| 2 + |Vx 2 | 2 ) = tr(cr(A)), 

where the second inequality follows from (14. 7j) . This is exactly (14. 10|) . 
To prove (|4.9|) . we will show 



^llmlLoo + jpllnihoo ^ c l- q / 4 , (4.12) 

for some constant cq = co(a), independent of A and L. We will subsequently adopt the convention 
that c is a constant, depending only on a, which can change from line to line. 

We first bound t\\. Let Q = (—1, l) 2 be the fundamental domain of the fast variable. Let d v Q 
and dhQ denote the vertical and horizaondal boundaries of Q respectively. Since Xi(?/l) 2/2 ) is odd in 
yi and even in y2, by symmetry we have xi = on d v Q, and d y2 x\ = on dhQ- Now if we consider 
the function xi + 2/1 > we have 

-A y (xi + y x ) + Av ■ V y {xi + yi) = 0, 

d 

|xi(y) < 1 on d v Q, and — (xi + yi) = on d h Q. 

Thus the Hopf Lemma implies X i + Vi does not attain it's maximum on dhQ, except possibly at 
corner points. So by the maximum principle X i + Vi attains its maximum on d v Q, and so 

llXlllioc < 1. 

Since X2 is bounded similarly, we immediately have 

I|Ti1 ] Il °° < cL- 1 < cL- a ' A . (4.13) 

The last step is to prove a bound on ||ri2|| i0 o- The crucial idea to bound r 12 is to split the 
right hand side of f)4.8j) into terms which are small in L p , and terms which can be absorbed by the 
convection term. To this end, write t\2 = rj + tpi + ^2 where rj, ipi are mean-zero, periodic solutions 
to 

2 

-A y r] + Av ■ V y r] = -2^2 d Vi Xi 

i=l 

-Aytpx + Av ■ Vyipi = A Y1 v i(xi + yi-^ sign(yj)^ - (viXi) 

-Ayip 2 + Av ■ Vyip 2 = -A Vj (yj - ^sign(yi) j . 

i=i ^ ' 

Before estimating each term individually, we pause momentarily to explain this decomposition 
of T\2- The equation for rj is of course natural. The equation for ipi stems from the well known 
behaviour of the corrector X l- We know from [6ll9lll9j that xi grows rapidly in a boundary layer of 
width 0(1/V~A) and decreases linearly in the cell interior. That is, for i = 1,2, 

1 
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away from the boundary layer. Further, by symmetry, Xi is odd in and even in the other variable. 
Thus, we expect the term Xi + Vi — \ sign(yi) to be away from zero only in the boundary layer (see 
Figure [3]), and hence should have a small L p norm! Now the equation for ^2 is chosen to balance 
the remaining terms, and thankfully the right hand side can be absorbed in the convection term. 

With this explanation, we proceed to estimate each function individually, starting with rj. Since 
V • v = 0, Lemma [27T1 guarantees 

IMIz,« < c II v xiIIl2. 

We remark that while Lemma 12.71 is stated for homogeneous Dirichlet boundary conditions, the 
proof in [4] goes through verbatim for periodic boundary conditions, provided, of course, we assume 
our solution is mean- zero. This justifies the application of Lemma 12.71 in this context. 
Since we know from [6] that ||Vxi|| L 2 = 0(A 1 ^' i ), we immediately obtain 

ML- < cA 1 ' 4 . (4.14) 

Our bound for tpx is similar in flavor. Let = Xi + Vi ~ \ sig n (yi)- Then for any p ^ 1, we know 
from [6] (see also |19} Theorem 1.2]) that 

Since —A(viXi) = 0(v~A), from Lemma 12.71 we have 

< cA (^\M LP + lte)IJ < c ^~^ ( 4 - 15 ) 

for any p > 1. 

Finally for ife, our aim is to absorb the right hand side into the drift. For i = 1,2, Let fa = fa a 
be defined by 



19 



and extended to be a C 1 , periodic function on M 2 in the natural way. Set 9 = ip2 + Yli=i(fi ~ (/»))> 
then is a periodic, mean-zero solution to 

2 

- Aj,0 + Av ■ V y = (M9i ~ A a /i) . 
i=i 

where ^ 

fl"i(l/) = d vJi ~ V% + 2 si S n (2/i) 
Since for any p ^ 1, we can explicitly compute 

\\gi\\ LP < cA~"k and ||A y /-|| LP ^ cA*~^ , 

by Lemma 12.71 we obtain 

ll^zIL- < l + ||%oa <cA 1_ £ (4.16) 

for any p > 1. 

Thus combining ([CH]) . ([USD and dSBD , we see ||ri 2 || iD o < cA 1-1 ^ 2 ^. Thus using (|1.14p and 
choosing p = when < a < 8/3, and p = oo for a ^ 8/3, we see 



8-3a 

IIT12|Il° 



L 2 

proving (|4.12|) . This completes the proof. □ 

Exit time from a square. 

Theorem 11.41 follows immediately from Proposition 11.41 

Proof of Theorem 11.41 The proof of Theorem 11.41 is now trivial. We simply inscribe a disk D = 
{\x\ ^ L/2} into the square D = [— L/2, L/2] 2 , and circumscribe a bigger disk D = {\x\ ^ L/y/2} 
around D. The corresponding exit times satisfy the inequality 

r(x) ^ t(x) ^ t(x), for all x G D. 

Using the bounds obtained from Proposition 11.51 applied to D and D, the inequality ()1.20j) follows. 

□ 

5 The eigenvalue in the homogenization regime 

5.1 The lower bound 

The lower bound for the eigenvalue stated in Theorem 11.21 follows . quickly from the upper bound 
on the expected exit time. 

Proof of the lower bound in Theorem 11.21 We claim that in general, we have the principal 
eigenvalue and expected exit time satisfy 

X^—^—. (5.1) 
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To see this, pick any e > 0, and suppose for contradiction that A < l/||r + e|| L00 . Then, 

-A(t + e) + Av ■ V(r + e) = 1 > ^ —r. (r + e) . 

||r + £||loo 

Also ^ 

— Ay? + • V<p = Ay? ^ ii |i — </j. 

Rescaling if necessary to ensure ||<y7||^oo ^ e, we see have ip ^ t + e in D. Thus Perron's method 
implies the existence of a function eft such that 

-A0 + AvVc/)= - —. i n 

||r + q| L oo 

= on dD, 

<P ^ ^ t in D 

This immediately implies 1/||t + e||^oo equals the principal eigenvalue A, which contradicts our 
assumption. Thus, for any e > 0, we must have A ^ 1/ 1 1 t - + ell^oo. Sending e — > 0, we obtain (|5.1D , 
Applying Theorem 11.41 concludes the proof. □ 

5.2 The upper bound 

In this section we prove the upper bound in Theorem II. 21 We will do this by using a multi-scale 
expansion of a sub-solution. As we have seen in the preceding sections, our multi-scale expansions 
are all forced to use a quadratic 'slow' profile, in order to avoid extra terms in the expansion. This 
makes the construction of the sub-solution slightly more difficult. As customary with homogenization 
problems, we rescale the problem so that the cell size goes to 0, and the domain is fixed. 

Lemma 5.1. Let h > 0, and ip be the solution of 

f -Aip + ALv(Lx) ■ VV = XB^ h in B 1 

[ V = ondB 1 , 1 ' ' 

where B r = {\x\ ^ r}, and xs is the characteristic function of the set S. Assume that A and L vary 
such that ()1.14p holds. Then there exists h > 0, and c = c(h) > such that 

c 

i/j(x) ^ -= for all x G Bi_ h . (5.3) 
V A 

provided A and L are sufficiently large. 

Lemma 15.11 immediately implies the desired upper bound. We present this argument below 
before delving into the technicalities of Lemma 15. II 

Proof of the upper bound in Theorem II. 2L Let <pi, fi± be the principal eigenfunction and the 
principal eigenvalue respectively for the rescaled problem 

-A(p% + ALv(Lx) ■ Vipi = fii(f\ in B\ 

if! = on dBi, (5.4) 
ipi > in B\. 
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Assume, for contradiction, /xi > where c is the constant in Lemma EH then 



V A 

— Aipi + ALv(Lx) ■ Vipi = ^ tpi. 

c 

Also, if ip is the function from Lemma 15. 11 then by the maximum principle, rp > in B\_h. Hence, 



'A 

-Atp + ALv{Lx) ■ Vip = Xb^ < ^r^- 

By the Hopf lemma, we know 5^ < on t?£?i, and so 991 can be rescaled to ensure ipi ^ V- Perron's 
method now implies that there exists a function <p that satisfies 



vA 

-Ac/) + ALv(Lx) • V(f> = —4> in Si 

= onaBi, ^' 5) 

> in B x , 



and, in addition, tp(x) ^ <^(x) </?i(x)- Therefore, fi\ = v A/c is the principal eigenvalue, which 
contradicts our assumption fj,\ > y/A/c. Hence fi\ ^ \J~Ajc. 

Now rescaling back so the cell size is 1, let A' and <p' be the principal eigenvalue and principal 
eigenfunction respectively of the problem (jl.lip on the ball of radius L/2. Since A' = 4/xi/L 2 , we 
have A' < 4\/A/(cL 2 ). Finally, let D be the square with side length L, and A, <p are the principal 
eigenvalue and eigenfunction respectively of the problem (II. lip on D. Then, since -B/,/2 C D, the 
principal eigenvalues must satisfy A ^ A', from which the theorem follows. □ 

It remains to prove Lemma 15. 11 

Proof of Lemma 15.11 Let 8 = t\ — ip, where t\ is the solution of (14. ip . the expected exit time 
problem from the unit disk. Now rescaling (jl.2ip to the ball of radius 1, (or directly using (|4,lip . 
which was what lead to (11.2ip ). we obtain 



< -^—=— 5.6 

VA 



Tl(x) ; — — - (1 — Ixl 2 

provided A and L are large enough and satisfy (|1.14p . Here C3 > is a fixed constant independent 
of A and L. Thus, (|5.3p will follow if we show that 

Ph-iB^^il-z') inf nix) (5.7) 

for some small e' > 0. Observe that (|5.6p implies that the right hand side of (15 .7p is Oihj \J~A). 
Therefore, to establish (|5.7p . it suffices to show that there exists constants ho > and c > such 
that for all h ^ ho, there exists Aq = Ao(h) an d Lo = Lq(K) such that 

Ml^b^) < 4f?l\ (5.8) 



provided A ^ Aq, L ^ Lo and (|1.14p holds. Above any power of h strictly larger than 1 will do; our 
construction below obtains h?l 2 , however, in reality one would expect the power to be h? . 
We will obtain (15. 8p by considering a Poisson problem on the annulus 

Al- 2h ,l = {1 - 2h ^ \x\ O}. 
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If we impose a large enough constant boundary condition on the inner boundary, the (inward) normal 
derivative will be negative on dB\_2h- Now, if we extend this function inward by a constant, we 
will have a super-solution giving the desired estimate for 6{x). We first state a lemma guaranteeing 
the sign of the normal derivative of an appropriate Poisson problem. 

Lemma 5.2. There exists Iiq and C2 > 0, such that for all h < ho, there exists Aq, Lq > such 
that the solution d\ of the PDE 

- A0i + ALv(Lx) • V0i = XAi A - h ™ M-2h,i 

6>i = on dB 1 , fo g\ 

1 = ^ 2 =/i 3 / 2 ondB^ 2h , 
VA 

satisfies 

— ^ on dBi_2h, 
or 

provided L ^ Lq, A Aq and (|1 . 14|) holds. Here ■Sp denotes the derivative with respect to the radial 
direction. Moreover, the function 0% attains its maximum on \x\ = 1 — 2h, and 9\(x) ?S C2/1 3 / 2 / \f~A 
for all x £ A\-2h,i- 

Now, postponing the proof of Lemma |5.2| we prove (|5.8|) . Choose h small, and A, L large, as 
guaranteed by Lemma 15.21 and define by 

h(x) when \x\ ^ 1 — 2h 

° {X) = < -^L/i 3 / 2 when |x| < 1 - 2h. 
I~A 1 1 



1 when \x\ ^ 1 — h 

when |z| < 1 — h k, \x\ ^ 1 — 2h 



where 0\ and C2 are as in LemmaEU Then £ C(#i) n C 2 (B^2h U ^i-2A,l), and 

(-A + ALv(Lx) ■ V)0» 
Further, when |x| = 1 — 2/t, 

* and -^<0 



<9r 9r 

where the second inequality follows from Lemma 15.21 Thus is a viscosity super solution to the 
PDE 

( -A0 + ALv(Lx)-V9 = XA 1 , 1 - h in B x 
\ = onaBi. 

By the comparison principle, we must have 0^0, which immediately proves (|5.8|) . From this (|5.7p 
follows, and using (|5.6p we obtain (|5.2p . concluding the proof. □ 



It remains to prove Lemma 15.21 Roughly speaking, if we choose the constant C2 sufficiently 
large, the function Q\ is nearly harmonic. The inhomogeneity of the boundary conditions dominates 
the right side of the equation. A "nearly harmonic" function should attain its maximum on the 
boundary, implying the conclusion of Lemma 15.21 

The reason we believe the constant C2/i 3//2 /V^4 is large enough, is because the homogenized exit 
time from the annulus is quadratic in the width of the annulus. Unfortunately, the slow profile 
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is not quadratic in Cartesian coordinates, and so the best we can do is obtain upper and lower 
bounds, which need not be sharp. We begin by showing that the expected exit time from an 
annulus of width h grows like /i 3//2 . While we certainly don't expect the exponent 3/2 to be sharp, 
any exponent strictly larger than 1 will suffice for our needs. 

def 

Lemma 5.3. Let A\-h,i be the annulus Ai-hi = B\\Bi_f l , andT ann be the solution of the Poisson 
problem 

f -Ar ann + ALv(Lx) ■ Vr ann = 1 in Ai- h ,i _ ^ 

\ T ann = on dAi-h,i 

Suppose L and A vary so that (|1.14p holds. Then there exists constants ho > and c > such that 
for all h < ho, there exists Aq = Ao(h) and Lq = Lq[K) such that 



provided L ^ Lq, A ^ Aq and ()1.14[) holds. 

Proof. The main idea behind the proof is that as A, L — > oo, we know that t\ tends to an explicit 
(homogenized) parabolic profile and is constant on dBi-^. Now if we subtract off a harmonic 
function with these boundary values, then we should get a super solution for r ann . Finally, we will 
show that a harmonic function with constant boundary values grows linearly near dB\, at the same 
rate as t\. Thus the above super solution will give an upper bound for r ann which is super dinear in 
the annulus width. 

We proceed to carry out the details. Let rj' be the solution of 



-An' + ALv(Lx) ■ V??' = in Ai- h ,i 

7]' = on dB\ 

' 2h ~ h2 f) R 

r?= 2lrMI)) OUdBl - h > 



(5.11) 



and define 



n - 77 + - 



where C3 is as in (15.6H . Then r ann satisfies 



-Af ann + ALv(Lx) ■ Vf ann = 1 in A\-h,l 

T ann ^ Oil dB ± 
Tann ^ On dB X -h- 

The first boundary condition follows because both 7/ and t\ are on dB\. The second follows 
from f|5.6[) and the boundary condition for 7/ . Thus, the maximum principle immediately implies 

that T an n ^ ^"aim- 

Since (|5.6p gives the asymptotics for t\, to conclude the proof we need a lower bound on 77' that 
is 'linear' in the radial direction near dB\. We separate this estimate as a lemma. 

Lemma 5.4. Let r] be the solution of 

-At? + ALv(Lx) • Vt? = in Ai- h ,i 

77 = on dB 1 (5.12) 
7/ = h on dBi-h, 
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Then there exists a constant c, independent of h, A and L, such that 

L -a/A\ 



7](x) ^ 1 



+ 



Vh 



(5.13) 



when L and A are sufficiently large. 

Returning to the proof of Lemma 15.31 we see that Lemma 15.41 gives 



rf(x) ^ 



h 



1 



2tr(a(A)) 

Now using the above and (jl.8p . it follows that 

cL" Q / 4 

Taring) < T ann (*) = T\(x) - T) (x) + 
1 



h 3 / 2 



+ 



Vh 



2ti(a(A)) 
1 

2tr(a(A)) 
1 



1 



(2 -/*)(! 



Vh 



(1 - |*|) (ft - (1 - |*|)) + c/* 3 / 2 ) + I J 



1 / cL~ Q / 4 



2\x(o{A)) 
obtaining Lemma 15.31 as desired. 

To complete the proof of Lemma 15.31 we need to prove Lemma 15.41 We do this next. 



□ 



Proof of Lemma 15.41 We will construct a sub-solution of equation (|5,12p in small rectangles over- 
lapping A\-h,i- For convenience, we now shift the origin to (—1,0), and consider new coordinates 

(xi,x' 2 ) = f (*i + 1,*2)- In these coordinates, let R be the rectangle of height 2/i 3//4 , width h and 
top left corner (po, h?/ A ), where po = 1 — (1 — /i 3 / 2 ) 1 / 2 (see Figure H|). 
We will construct a function ry su b such that 



-Ar] snh + ALv(Lx) ■ Vr? sub = in Ai- h ,i, 

%ub ^ h on the right boundary of R, 

Vsub ^0 on the other three boundaries of R, 



(5.14) 



and Tfeub satisfies the linear growth condition 



%ub(4,0) ^x' 1 -c[h^ 2 + 



L -a/4 

Vh 



when x' x £ [po, p + h], 



(5.15) 



for some constant c independent of L, A and h. 

Before proving that the function 7] su ^, exists, we remark that by the maximum principle, we 
Vsub ^ V on R n Ai-hh- Moreover, as po = 0{h 3 / 2 ), the estimate (pU3l) can be extended to 
x'i £ [0, po] as well, possibly by increasing the constant c. This proves Lemma 15.41 when * is on 
the negative xi-axis. Now, if ()5. 15j) is still valid when the coordinate frame is rotated, our proof of 
Lemma 15.41 will be complete! 
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Figure 4: Domain for 77 su b- 



We will first prove that a function 7? su b satisfying (|5.14p and (|5.15p exists. We will do this by a 
multi-scale expansion. Let 

1 1 CqL 

Vsub(x') = rjo(x') + -771 (a/, y) + j^miv) ^/%~' 

where y = Lx' is the fast variable, tjq is given by 

t i i \ i , (( x 'i) 2 ~ ( x 2) 2 ) 

and Co, and ci are constants, each independent of L, A and /i, to be chosen later. As before, r\\ is 

2 



i=i 

and r/2 is the mean 0, periodic solution to 

2 



A y r) 2 + Av(y) ■ V y r] 2 = ^ [%d yi Xj - A (viXj - (viXj))} d x >d x >m, (5.16) 



3 

1 = 1 



where Xi are the solutions to (|4.7p . 

Again, the crucial fact here is that since 770 is quadratic, the second derivatives are constant and 
772 becomes independent of the slow variable x' . Using this, a direct computation shows that 

- Ar7 sub + ALv ■ Vr] suh = -Ar/o + A ^(^j)^^?jo. (5.17) 

Note that by symmetry, (vixi) = {V2X2), and (7J1X2) = (^2X1) = 0- Further, by our choice of 770, we 
have Ar7o = 0. Hence, the previous equation reduces to 

-Ar7 sub + ALv • V77 sub = 0, 
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as required by the first equation in (|5.14j) . 

Next, we show that if c\ is appropriately chosen, we can arrange the boundary conditions claimed 
in (|5.14j) for rjo- Notice that po = 0(/i 3//2 ), and on the top and bottom boundary we have x' 2 = ±/i 3 / 4 
and x[ G (po, Po + h). Thus 

4 and K^O^ 3 / 2 ). 

So choosing c\ large enough, we can ensure 770 ^ on the top and bottom of R. 

On the left of R, we have x[ = po and \x2\ ^ /i 3//4 . So (x'^/Vh ^ 0(/i 5 / 2 ) = o(po), and choosing 
ci large we can again ensure 770 ^ on the left of R. Finally, on the right of R, we have x\ = po + h 
and \x' 2 \ ^ /i 3//4 , and we immediately see that for c\ large enough, we have 770 ^ /i on the right of R. 
Thus 7/0 satisfies the boundary conditions in (15.14p . 

To see that r7 su b also satisfies the boundary conditions in (|5.14|) . we need to bound the correctors 
appropriately. Exactly as in the proof of Lemma 14.11 we obtain 



cL~ a / 4 

< - vr- (5 - 18> 



where c > is independent of L, A and h. We remark that the extra \j\f~h factor arises because 
derivatives of 770 are of the order 1 j\ph and they appear as multiplicative factors in the expressions 
for 771 and 7/2. 

Consequently, if cq is chosen to be larger than c, we have 

1 1 c L- Q / 4 

I" 1 + ~ *■ °' (5 ' 19) 

Since 770 already satisfies the boundary conditions in (|5.14p . this immediately implies that 7/ su b 
must also satisfy these boundary conditions. Finally, since 770 certainly satisfies (|5.15p . it follows 
from (|5.19p that 77 su b also satisfies (|5.15p . This proves the existence of 7/ S ub- 

Now, as remarked earlier, the only thing remaining to complete the proof of the Lemma is to 
verify that if the rectangle R, and the coordinate frame are both rotated arbitrarily about the center 
of the annulus, then there still exists a function 7/ su b satisfying (|5.14p and (I5.15P in new coordinates. 
This, however, is immediate. The new coordinates can be expressed in terms of the old coordinates 
as a linear function. Consequently, our initial profile for 770 will still be a quadratic function of the 
new coordinates. Of course, by the rotational invariance of the Laplacian, it will also be harmonic, 
and the remainder of the proof goes through nearly verbatim. The only modification is that after 
the rotation the mixed derivative d x i d x i rjQ no longer vanishes, and the terms involving V1X2 and 
V2X1 do appear in (|5.16p and (|5.17p . However, they can be treated in an identical fashion, as in 
the proof of Lemma 14.11 using the precise asymptotics for xi an d X2 from |19j . This concludes the 
proof. □ 



Finally, we are ready for the proof of Lemma 15.21 
Proof of Lemma 15.21 For a given C2, let ij' be the solution of 

-At/' + ALv(Lx) • Vt/ = in Ai- 2 h,i 



77 = on dB\ 

'A 
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Then, we have 

01 - rf = on dAi-2h,i, 

and 

-A6>i + ALt>(Lz) • V(9i ^ 1 in ^ti_ 2 / l ,i- 

Consequently 6\ — rf ^ T ann , where r ann is the solution of (j5.10j) on the annulus A\-2h,i- Thus 
applying Lemma |5\3| we see 

9i(x) < — L(2/i) 3/2 + f/(a?) (5.20) 
V A 



The function r/ decreases at most linearly with \x\. This can be seen immediately from an 
asymptotic expansion for a super solution. Indeed, starting with 

m{x) = ^/I {l - Xl) 
and choosing r]i and 772 as in the proof of Lemma 15.41 we immediately see that 

rf( Xu Q) ^ ^(i_ Xl ) +c ^- 1 /2 L -«/4. 

We remark again that the extra \fhA~ 1 / 2 factor arises from the gradient of 770 • Now by rotating the 
initial profile 770 appropriately, we obtain the linear decrease 

rf( x ) < ^2(1 - Id) + cVhA- l ' 2 L- a /\ (5.21) 

as claimed. 

We claim that (|5.20|) and (15.21j) quickly conclude the proof. To see this, note first that equa- 
tion (|5,20p and (|5.2ip immediately give 

01 fx) ^ ^=(2hfl 2 + -^=/i 3 / 2 + cVhA-^L- / 4 whenever x E A^ h 1. (5.22) 
y/A 2VA ' 



However, since 

-A0i + ALv(Lx) • V6»i = 

on Ai—2h i-hi the maximum principle implies that 9\ can not attain it's maximum in the interior of 
the annulus Ax-2h,i-h- Consequently (I5.22j) must hold on the interior of the entire annulus A\-2h,i- 
Now if we choose C2 large enough so that 2 3 / 2 ci < ^, and then choose L, A large enough so that 
cA- l / 2 L- a / A < f , we see that 81 is forced to attain it's maximum on the inner boundary dBi_ 2 h, 
and the Lemma follows immediately. □ 
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